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Abstract 

We derive the stochastic Schrodinger equation for the limit of contin- 
uous weak measurement where the observables monitored are canonical 
position and momentum. To this end we extend an argument due to 
Smolianov and Truman from the von Neumann model of indirect mea- 
surement of position to the Arthurs and Kelly model for simultaneous 
measurement of position and momentum. We only require unbiasedness 
of the detector states and an integrability condition sufficient to ensure 
a central limit effect. Despite taking a weak interaction, as opposed to 
weak measurement limit, the resulting stochastic wave equation is of the 
same form as that derived in a recent paper by Scott and Milburn for the 
specific case of joint Gaussian states. 



1 Introduction 

The theory of continuous measurements of a quantum system, based on the 
Ludwig formalism, was originally presented by Barchielli et al. 1 . Since then 
stochastic Schrodinger equations describing the dynamical evolution of the state 
of a system conditional on observations of a monitored set of observables iXj j 

have been developed by several authors [211300 E! The generally accepted 
form is an adapted stochastic differential equation for a vector state valued 
process |\T/) of the type (see, for instance [7]) 



m dt- > x/2k7 I X,- - ( X,- > ) I*) dB t 



,00 

(1) 
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where (Xj) = y&lXj^j are current observations (thus making the equation 
non- linear) and {B^} is a multidimensional Wiener process. The constants k 3 
are positive and H is the Hamiltonian governing the background evolution. 
Note that the stochastic wave function is normalized, since, by the Ito rule 



dB[ 3) dB { t k) = Sjkdt, 



d ||* |p = (d#|#) + (#|d#) + (d*|dW) = 0. 
If Y is an observable and if we set ( Y ) = ( \&| Y^f) then 



d ^ = (C (y)} dt - cov (y, X 3 ) dB 



GO 



(2) 



where £ is the (self-dual) Lindblad generator 
1 



C Y = 



Y, H 



Xj,Y 



Xj + Xj 



Y,Xj 



and we have the covariancc 



(y,x) = 1 -(xy + yx)-(x)(y) 



Stochastic Schrodinger equations can alternatively be derived [S] from quan- 
tum stochastic calculus 0. Recently, a probabilistic derivation based on limit 
theorems was presented by Smolianov and Truman ^Uj which considered re- 
peated position measurements according to the von Neumann model. In this 
paper, we recap on the calculations in |1(J| and correct some numerical errors. 
We then generalize their argument to consider the simultaneous measurement 
of position and momentum according to the Arthurs and Kelly model this 
is non-trivial as we have to work with quantum rather than classical probability 
distributions. This model has been examined by Scott and Milburn ^21 with 
explicit calculations made for Gaussian states of the measurement apparatus. 
They obtain the crucial result that the resulting stochastic wave equation for 
continuous limit of weak measurements will be of the form (|TJ with X\ = q and 
X2 = p. Our approach uses a different limit: a weak interactions limit. (Our 
treatment of the Smolianov and Truman shows that the limit can be interpreted 
as either as finite interactions with weak measurements (as they consider) or as 
finite measurements with weak interactions: actually both limits yield the same 
result.) We consider a general classes of states for the apparatus where the only 
restriction is unbiasedness for the pointer positions and momenta and a sim- 
ple integrability condition, that the constants k given by l|13|l below are finite, 
which ensures a central limit effect. The resulting stochastic wave equation is 
nevertheless of the same form as that obtained by Scott and Milburn. 
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1.1 Measurement Conditioned Wave Functions 

We consider a fixed Hilbert space f)s = L 2 (M n ) which describes our systems and 
consider indirect measurements made on a second system (the apparatus) having 
Hilbert state space and t)A = L 2 (R m ). The joint space is l)s^^)A- A vector state 
<f> can be considered as a function (f> = <p (q; q') where q and q' = (q[, . . . ,q' m ) are 
appropriate coordinate variables for the system and apparatus respectively. Let 
(q[, . . . ,q' m ) denote the corresponding position observables for the apparatus: 
they form a maximal set of commuting, self-adjoint operators on t)^. 

Suppose that the system and apparatus are initially prepared independently, 
so that the joint state is the vector O = <p s ® <f> A . Subsequently, the two un- 
dergo an interaction described by a unitary operator V. The state immediately 
prior to measurement is then <f> — V (f> . We then measure the observables 
(q[, . . . , q' m ). According to the rules of quantum mechanics, we obtain a family 
of classical random variables Q' = (Q[, . . . , Q' m ) describing the measured values 
Here Qj is the random variable describing the measured value of the observable 
q'j. distributed according to 



E 



exp 



3 = 1 




j 

Jr> 



p(q')exp<jz^^ g ^d"V, (3) 

3=1 



where p(q') = /|(q;q'|0)| 2 <Pq. 
We then consider the function 



V(q|q') := 



<q;q'l0> 



and introduce the stochastic wave-function VP defined by 



(4) 



<q|*)=V(q|Q') (5) 

If we do not wish to make the q-dependence explicit, we shall write \& = ^q/. 
The stochastic wave-function interpretation comes about as the mathematical 
re- interpretation of the wave- function as the \) 5- valued random variable depen- 
dent on the measured random variables Q' = (Q[, . . . , Q' m ). 
Let X is an observable of the system, only. We set X (Q') = 

ftp* (q|Q') (x^j (q|Q') d n q. Then the average value of the random variable 
X (Q') will be 

E[l(Q')]=|l(q')p(qW (6) 
and this is clearly ($\X = f <f>* (q;q') (jt<j?) (q;q') d n qd m q'. 
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1.2 Repeated Measurements 



Next suppose that several measurements are made at regular intervals of time 
r. Let A (j) denote the apparatus employed for the j-th measurement and let 
4>Atj\ be its state prior to interaction with the system. (We may either consider 
an assembly of apparatuses A (1) , A (2) , A (3) , • • • or a single apparatus whose 
state is forced to be <j>A(j) somehow just before the j-th measurement for each 
j.) Between measurements, the system is allowed to evolve according to the 
Hamiltonian H on f)g. 

If we denote the conditioned wave-function immediately after the n-th mea- 
surement by ip n = ip n (x\Q' 1 ; - ■ ■ , Q^), then we have the iterative relation 

1& n (q|qi;--.;q;) = -^= U <£>\ e T&/ih V n V>„-i ® • (7) 



2 Measurement of Position 

We begin with the simplest problem of monitoring a single observable. We take 
this to be a position coordinate q. 

As system-apparatus interaction, we take a linear coupling of the system 
position with the apparatus momentum p' . The interaction is then given by 
the unitary V = exp {fiqp' / ih} . We then measure the position observable, 
g', for the apparatus. This is essentially the von Neumann model for indirect 
measurements. The apparatus state prior to measurement is taken to have a 
real symmetric wave-function and so its mean position and moment are zero. 

The calculations in this section follow the same arguments as in Smolyanov 
and Truman |10j . 



2.1 Von Neumann's Model 

The interaction between the system and apparatus is then given by 

<Ps (?) 4>a (</) <Ps (?) 4>A W ~ W) ■ (8) 
The probability density for the apparatus position after interaction is then 

P(q') = J \<P S {q)\ 2 \cj> A {q' ~m)\ 2 dq (9) 

which is a convolution of two probability densities. 
We choose (j> A to have the form 

^<«'> = ^p) < 10 > 

where x is real-valued and normalized so that x (y 2 ) dy = 1 = x (j/ 2 ) V 2 dy- 
In particular, Y = —Q' will be a mean-zero, unit-variance random variable. We 
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then have the decomposition 



Q' = aY + ft (Z + (q)) 

where Z = Qs — (q) is centered. 

Assuming that \ is analytic, we get an expansion of the form 



(11) 



1 



4>a (Q'-w)-^x(^ 2 ) jl + 

where s = a~ 2 \Q' — pq\ 2 — Y 2 . 



X '(Y 2 ) l X "(Y 2 ) 



-e + 



X (Y 2 ) 2! X (Y 2 ) 



Likewise the factor p (Q r ) x ^ 2 can be expanded. First of all we observe that 
J \4> s (q)\ 2 edq = 2^-ZY + (£) * (Z 2 + a 2 q ) 
j\<i> s {q)\ 2 e 2 dq = 4^)\z 2 + „ 2 q )Y 2 + 0^j 



where we set a 2 := (q 2 ) — (q) 2 (the variance of the observable q for state cj> s ) 
and treat (f ) as small parameter. 



p(Q') = \<f>s(q)r-oX 



2 1 2 / \q' - M?r 



dq 



fu I .,2 X 2 (y 2 ) f. , 2 X '(Y 2 ) 1 



X'(Y 2 ) 2 , x"(Y 2 ) 



X (Y 2 ) 2 X{Y 2 ) 



X 2 (Y 2 ) \ ^(Y 2 ) p 



+4 



X'(Y 2 ) | 2 X'(Y 2 ) 2 | 2 x"(Y 2 ) 



X(Y 2 ) ' " X (Y 2 ) 2 ' " X(Y 2 ) 
We obtain to lowest orders 



(Z 2 + a 2 q ) Y 2 + o((t) 



i>{q\Q') 



')={l-2^(g-<g»r+ 



+ 



X i o X- \r2 



X X 



+ 2^Y 2 I (g - (q)Y - I — + 2— r 2 + 2 



^ 2 <T 2 



X 



+2Z [2(^1 r 2 - _ 2 ^r 2 1 ( g - (q)) 



X X 



+0 ((a) 3 )r 5(<z) - (12) 



where x, x' and x" are evaluated at Y 2 . Note that the terms involving Z 2 cancel 
to this order. 
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2.2 Continuous Measurement Limit 



We now consider sequential measurements with the successive apparatus states 
prior to interaction being copies of (j> A chosen as above (I1U|) . 
Adopting the scaling (f ) = y/r, we get 



ih 



j» + 2^) iq - , syr - 

Xn Xn 



TXr, 



+ 2^K? + 2 



Xn 



A, n 



Xn Xn / 



> An t/"2 



+ o(r 1 / 2 )U n _ 1 . 



where x n = X (Y 2 ), etc. 

In the limit t — > we expect from the law of large numbers that 



^ fx'- y" 



Xj 

[t/r] 



-Kt; 



3 = 1 



for suitable /. Here k — — + 2^y 2 ^j x 2 dy, after some algebra and an 

integration by parts we have that 



(x'x + 2xxV ) dy 



(x'y) dy>0. 



(13) 



Inspecting the coefficient of a v we are lead to consider 

2 



[t/r] 

lim t > 

T->Q ^ 



X 



-1 + 2 *Ly/ + 2 

Xi Xj 



II 



>7 



6>i; 



where 6* = /_ (xx' + ^x' 2 y 2 + %x"xy 2 ) dy and comparison with the above 
integration by parts shows 9 = 0. 

Likewise, from the central limit theorem, we expect that 



lim 1 V X ' ( ^V.- fE 



where B t is a standard Wiener process. 



We therefore are lead to the stochastic differential equation 
ih 



\<M) = i--H-K(q-(q)) 2 \\y) dt- ^2 K (q - (q)) |tf ) dB t . CI 4) 
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Remark 1 The limit £ — > can occur in two distinct ways. We can consider 
a weak coupling limit fj, — ► where the measurement process is fixed (<r = 1). 
Alternatively we can consider a weak measurement limit a — > oo where the 
interaction is fixed (fj, = 1). 



Remark 2 27ie equation \14\j can clearly be generalized to QJ) provide that 
consider indirect measurement of commuting system observables. 



we 



Remark 3 The result \1$ corrects some numerical errors made in JflClj/ with 
regard to the coefficients. 

3 Measurement of Position &; Momentum 

We take the system to have canonical position and momentum observables q 
and p, respectively. Our apparatus consist of two distinct components, A' and 
A", which have canonical observables q',ft' and q",p", respectively. Our aim is 
to couple the system to the apparatus by interaction and make inferences about 
the system based on simultaneous measurements of the position of A' (that is, 
q') and the momentum of A" (that is, ft"). 

Initially we take the system to be prepared in state <p s while the apparatus is 
prepared in a state 4> A £ §a = §A' ® f)A"- We assume that the two components 
of the apparatus are not entangled - that is, 4> A (q',p") = 4>a> (sO ^a" (p")- 

The interaction between the system and apparatus is described by the uni- 
tary operator 

V = cxp {{nft'q - vq"p) /ih} (15) 
and we dehne an operator V (q',p") on t)s by 

(V (q\p") <t> S ) (q) = («, q',p"\V <j> s <j> A ) ■ (16) 

We see that 

V (q',p") = I ^f^f ( q >\ p >) ( p »\ q ») W (-uq", -rf) {p'\q') (q"\p") (q',p" 
= I (q' - M,p" - vp\<$> A ) W (q,p) 

and so we obtain V (q' ,p") as a Weyl-quantized operator (see appendix) 

V(q',p") = [{q'-pq,p"-vp\<t>A)]w ey i- (17) 

We next denote the random variables obtained by observing q' and ft" by Q 1 
and P" respectively. As the observables commute, we are able to assign a joint 
probability to Q' , P" once a density matrix is prescribed. 
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A stochastic wave-function dependent on the observed position and momen- 
tum variables of the apparatus is 



* = 



V{Q',P") 



where p(q',p")=J (q,q',p"\V(t> 



dq = 



V(q',p")<l>s 



(18) 



Under the action of the unitary V we have the relations 



VU'V = q' + M - \iiuq" 



V ] p" V =p" +vp- -fJ,vp'. (19) 



(The action of V can be understood as that due to the evolution in unit time 
governed by Hamiltonian [ip ' q — uq"p. This Hamiltonian generates a linear set 
of equation for the canonical variables which is readily soluble. The variables p' 
and q" are evidently invariants.) 

The joint probability distribution of Q' and P" is determined from the char- 
acteristic function 



,i(aQ'+l3P") 



(aq'+fjp")^ 



where <fi = V {<p s ® (j) A ) is the post-interaction state. Using the lemma, we can 
write the characteristic function as 

4> s \e l{a ^ +f}vi,) 4>s) x (4>A'\^ aS '~^^')(f> A ,) x (^«|e i ( w "-*°'*"*'')0 j4 „\ 

The result is that we have the following decompositions into sums of independent 
random variables: 



1 



1 



Specifically, 



' = Q'o + mQs - /^Qo'; p" = Pq + vP s - -ij-vPq 



Pr[q<Q s <q + dq] = \<j> s (q)\ 2 dq 



(20) 



\<f> A , {q')\ Z dq' 
|2 



Pr[g' <Q' <q' + dq'} 
Pv[q"<Q^<q" + dq"] = (</"k r }| J ,/,/' 
while a similar set of laws hold for the P's. 



Remark 4 Some caution is needed here: a joint distribution for Q' , P" exists 
as they correspond to commuting observables, however this is not the case for 
any of the pairs Qs, Ps or Q' Q , Pq or Qq, Pq . 

Remark 5 The parameters /x, v are free. In JT^j they are chosen as v = so 
that symplectic area is preserved: the single parameter fi is referred to therein 
as a squeezing parameter . 
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3.1 Perturbative Expansions 

We take the initial states of the apparatus to be 

<t>A'W)=x{W\ 2 ), 0a»(p") = a(|p"| 2 ) 

Here we shall adopt the weak coupling scheme and fix a = 1 while fi and v are 
the small parameters. Expanding to lowest orders, we see 



(q - fj,q,p" - vp\4> A ) = x (V - w| 2 ) A {\p" - vp\ 

where \ — X (l^'| 2 ) an< ^ A = A (jp"| 2 ^ , etc., and 

e q = -2[iq'q + fi 2 q 2 , e p = -2up"p + v 2 p 2 . 
Taking Weyl quantization yields 

V (q',p") = XA (l - 2^q'q - 1v^-p"p 
V + 2^ 9 ' 2 ) q 2 + v 2 + 2^>) p 2 + 2 A uX£ 9 V' (# + #) + -} 
and the joint probability density for the observations is 

p(q , ,p") = {<t>s\V(q',p"f<t>s) 

= X 2 A 2 (l-4 M ^' (q)-4^p" (p) 
{ X A 

<«f> + , 2 



+2^ 



X' , „X" , 2 \ , o / xY X 



+ 2^ +2, 

x x / V x 



A' A" „ 2 \ / A'p" N - 

A +2 T p +2 



(f) 



x' A' I 

+6/ZI/ -q'p" (qp + pq) H } 

X A J 

We therefore obtain the approximate form 

V(q',p") 1 9 X' 9 A' 

= 1 - 2^— <? (g - (qr)) - 2^— p (p - (p» 



VW/) X A 
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+6 



A 



X'A' , 



(P) -4 



Ay 



p(p) 



+6/iz/— ( 2 (SP + W) - 3 (<ZP + M) + 6 (g) (p) - 4g (p) - 4p (q)) + 

X A 



Our objective is to study the random propagator 



v(q',p") 



where now the 



classical random variables Q' and P" have p as their joint probability density 
function. As we have seen, the exact distributions of Q' and P" are non-trivially 
related to the original states of the system and apparatus. However, we are only 
interested in the lowest order dependence in terms of parameters fi and v. The 
decompositions <|20|) may be rewritten as 

Q' = Q'o + A* («) + liF, P" = Po+" (p) + vG (21) 
where Q and Pq are independent variables with the pre-interaction densities 

X (l?o| 2 ) an d A^|p '| 2 ^ respectively, and F and G are further independent 
variables with means of order v and /i, respectively. To lowest order, we obtain 

X' 



'(\Qf) 
' Q' 



(iQf) 
' (>"1 2 ) 

A(|P'f) 



A 



P" = 



Xo 



A " 



Xo 
Xo 

K 

A n 



2 Xo _ Xo 



Xo 



Xo, 



2 ^--(^ 
Uo Uo, 



Q'o 2 



p" 



(v2> + ^)- 
«P>+G) 



where x = X ( IQol' 
to 

V(Q',P 



A = A (\Pq\ ), etc. Making these replacements leads 



= i _ 2ai a4q (( j _ ($)) _ 2v^P{ (p (p» 
y/p{Q',P") Xo A 
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2^g? ) f 



Xo Xo 

2 



Xo Xo 



<<f>+6 



Wo 



V Xo 

a; 



+ 2 f Xo _ (XoY ] g'2 
Xo VXo VXo/ J ° 



\ Xo J J V Xo 



.Aft 



A ' 2 A P ° ]P 



, A' P 



(p 2 ) + 6 



A' P" x 2 



A 
Xo A 



P (p) - 2 



A lAo \A 



pn 
r o 



A y y x " 7 ■ " v A 

((p)+G)(p-(p)) 



Xo A o 

4 Repeated Measurements 

Next suppose that simultaneous measurements of position and momentum are 
made at regular intervals of time r. Again, we assume that the pre-interaction 
states 4>A{j) are an copies of the fixed state <f> A . 

If we denote the conditioned wave-function immediately after the n-th mea- 
surement by tp n = 2p n (x\Qi, P{'; • • • , Q' n ,P%), then we have the relation 

,/, - p rH/ih V (Q'n, P n) , 



Setting \i = v = y/r, we can use (1221) to obtain 

(23) 

To obtain Q from (J23) we made the replacements ^+2^Q 2 with -2 (j&Q'o) ^ 

and ^ +2 A ^-Pq 2 with —2 f^fo') which is consistent with (fTTSf) . In particular, 
the coefficients of (g 2 ) and (p 2 ) disappear. The terms involving F and G are 
dropped since they are O (t 1 / 2 ) and so will give negligible contribution in a law 
of large numbers limit. Finally the cross terms vanish since 



T ^ 0r ^Vxo 7,V A « 



ii 



on account of the fact that / \ (y 2 ) x' (y 2 ) ydy = = / A (x 2 ) A' (x 2 ) xdx. 



Comparison with position-only situation shows that we are lead to the 
stochastic differential equation 

|d*) = l [ ^H-K q (q-(q)) 2 -K p (p-(p)) 2 ^mdt 

?(«) _ FTT (A — iaW iiTi\ ^r(p) 



where 



-y/2 Kq {q - (q)) |*) dBl"' - ^/2k p (p - (p)) |*) d£ 



«, = 2^ (x (y 2 ) 2/) 2 dy, n p = 2 J (A' (a; 2 ) xf dx 



and we have the following limits in mean-square sense to independent Wiener 
processes 



B 



(?) 



[t/r] 



lim 

V K i T *=i ^ Xo 



? (p) 




hm 

1 3 = 1 



[t/r] 

E 



4.1 Stochastic Dynamics 

The associated Lindblad generator for this monitored dynamics is 

C(Y)=±[Y,H] + Kq {[q,Y] q + q[Y,q]}+K P {[p,Y]p + p[Y,p\}. 



We note that C Y 



^Y, occurs for the special case of observables of 

the type Y = aq + (3p + jqp. This means that when H = -^p 2 + $ (q) , the 
averages of the canonical observables evolve in a non-random way according to 
the Ehrenfest theorem for closed systems: 

£»- = »• !<#> = -<*'«»■ 

However, we do not have the derivational property C (XY^J = XL (y) + 
C[X)Y and so the dynamics is dissipative. In particular, 



C 



(H)=n q - + Kp hH"(q) 



and so energy is not conserved. 



4.2 Conclusion 

We have shown that the central limit effect allows us to derive a stochastic 
Schrodinger equation in a very general setting. Essentially we need only be 
in the domain of attraction for Gaussian statistics. (It is possible that more 
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general results hold for stable laws - leading to stochastic Schrodinger equa- 
tions driven by Levy processes - however, this clearly would be indicative of an 
imperfection in the measurement apparatus only.) The striking result is that 
the non-commuting observables can be measured simultaneously with negligi- 
ble interference. This backs up the general phenomenological approach used in 
areas of quantum control and filtering. In their article, Scott & Milburn [T2] 
investigate several models, including classically chaotic ones, and show that the 
system is localized with the monitored trajectory in phase space correspond- 
ing to a quantum average plus noise. Whereas there are still many interesting 
open questions regarding, for instance, the semi-classical limit, the present re- 
sult does show claim that the analysis of Scott & Milburn approach is generic 
for continuously monitored phase variables. 



5 Appendix 

Let (q,p) £ Y where Y = M 2 is phase space. For a given function f = f (q,p) 
on phase space, the association of an operator / (q,p) is ambiguous due to 
the problem of operator ordering. We shall adopted the Weyl quantization 
convention. 

The Weyl operator at phase point (q,p) is defined to be W (q,p) = exp (qp — pq) 
From the CCR and an application of the Baker-Hausdorff-Campbell theorem, 
we find that 

W(q,p)W(q',p') =e( qp '~ pq ')/ 2lh W(q + q>,p + p'). (24) 

Let f = f (q, p) be absolutely intcgrablc on Y and define its phase space 
Fourier (Weyl) transform to be 

/(9 ' P):= / r ^ e( " )M /feP")- (25) 
The Weyl quantization of / is then defined to be the operator 

[/(<7,j5)] W eyl : = J r ^W(q,p)f(q,p) 

= I ^^e^-^f(q,p)W(q,p) (26) 

and we refer to the map / (q,p) i— > [/ (<Z,p)] Weyl as Weyl quantization. 

We have for instance [exp (t (qp — pq) /ih)] Wcyl = W (tq,tp) and expanding 
in powers of t, we obtain [(qa + pf3) n ] Weyl = (qa+p(3) n . Further expansion 
in terms of a, /3 show that polynomials will be mapped to the symmetrically 
(Weyl) ordered form.. For instance, [q 2 p] Woyl = | (q 2 p + qpq+pq 2 ), etc. 
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